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We study the Hubbard-Holstein model using slave-boson mean-field and a variational Lang-Firsov 
transformation. We identify weak and strong e-ph coupling solutions, whose stability depends both 
on the bare e-ph coupling and on the correlation strength. At mean field level the evolution from 
weak to strong electron-phonon coupling occurs via a first-order polaronic transition if the adiabatic 
parameter is below a critical value. In the strongly correlated regime and in the adiabatic limit, the 
region in which the weak-coupling solution is stable is sizeably enlarged with respect to the weakly 
correlated system and the Mott metal-insulator transition is found to be robust with respect to e-ph 
interaction. 
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I. INTRODUCTION 

Experimental indications 1 of noticeable electron- 
phonon (e-ph) effects in the superconducting cuprates 
have revived the discussion of the interplay of e-ph and 
electron-electron (e-e) interactions. This issue has ob- 
viously a broader relevance, as both interaction mech- 
anisms are important also in materials like fullercncs 2 
and manganites*. The wealth of experiments available 
for the cuprates underscores quite a puzzling scenario 
where the fingerprints of e-ph interaction are strikingly 
different than in weakly correlated materials. In fact, 
while some quantities seem to be affected by phonons 
even more than in ordinary metals, particularly in the 
underdoped regime (photoemission spectra*, penetration 
depth and superfluid density**), others are almost in- 
sensitive to phonons, like resistivity. Such a situation 
is ultimately related to the presence of several energy 
scales in strongly correlated materials, which makes it 
compelling the use of nonperturbative approaches. Sev- 
eral nonperturbative methods have been indeed applied 
to different models with both e-e and e-ph interactions. 
Among those, we can count fully numerical approaches 
such as various kinds of Quantum Monte Carloi***, ex- 
act diagonalizatiori*** and Density-Matrix Renormaliza- 
tion Group**, the semi-analytical dynamical mean-field 
theory (DMFT)******, and approximate analytical re- 
sults based on large-N expansions****. Each of this ap- 
proaches has its own advantages, and typical range of 
validity. An analysis of the recent literature seems to 
uncover an unclear scenario, since some calculations con- 
clude an enhanced role of e-ph interaction in the presence 
of correlations)*!* and others suggest exactly an opposite 
situation. 14 An important part of those discrepancies can 



be attributed to the different physical regimes in which 
the investigations have been performed (e.g., antiferro- 
magnetic or paramagnetic phase), but another part is 
probably due to the numerical character of the previous 
analyses which does not elucidate enough the underly- 
ing physical mechanisms. In the present paper, we try 
to build a more analytical understanding for the corre- 
lated metal at and near half-filling and compare with the 
results obtained using DMFT in the same regime near 
the Mott transition. 14 DMFT is the quantum version of 
mean-field theories, that freezes spatial fluctuations but 
fully retains local dynamics. The method has been re- 
cently used to study the effects of e-ph interaction close 
to the Mott metal-insulator transition in the Hubbard- 
Holstein modelii*** In particular in Ref. 0, it has been 
shown that, close to the Mott transition, e-ph interac- 
tion manifests in a different way in low- and high-energy 
properties. While the high-energy Hubbard bands are 
affected by e-ph coupling, the low-energy quasiparticle 
physics is to some extent "protected" against the phonon 
effects by e-e correlation. The main effect of e-ph inter- 
action is indeed the reduction of the repulsive U leading 
to an effective interaction which can be parameterized in 
a simple way. Since DMFT still requires numerical calcu- 
lations which sometimes make it difficult to disentangle 
different effects, it is instructive to compare approximate 
analytical methods with DMFT in order to gain insight 
into the physical picture arising from the latter method. 

To this purpose we present here a mean-field picture of 
the Hubbard-Holstein model, following a recent sugges- 
tion by Perroni et alH. From one side we exploit the abil- 
ity of the slave-boson mean-field approach 18 to capture 
the main features of the Hubbard model in the strongly 
correlated regime. To treat the phonon degrees of free- 
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dom we use a variational Lang-Firsov transformation 19 , 
that has been introduced to describe the effects of c- 
ph interaction both in the weak- and strong-coupling 
regimes 20 . To this end one introduces a "rotation" pa- 
rameter /— , defined below, which measures the degree 
of polaron formation. 

Slave-boson approaches have been recently 
performed 21,22 , in a similar spirit of the present 
paper, to gain insight on the Quantum Monte Carlo 
results of Ref. (i. proposing that, in a range of corre- 
lation values, the small momentum e-ph vertex can be 
enhanced by increasing the repulsion (yet being smaller 
than the bare one). The studies of Refs. I2ll22l suggest 
that the enhancement of the vertex could be related to a 
finite-temperature signature of a zero-temperature phase 
separation. Our approach here follows the spirit of those 
calculations (namely the use of simplified approaches to 
gain insight into more involved calculations) but focuses 
on a slightly different subject, such as the "protection" 
of low-energy degrees of freedom from e-ph effects. 

The paper is organized as follows: in Section [n] we 
introduce the model and derive the general mean-field 
equations for the paramagnetic homogeneous phase. In 
Section IIIII after a brief summary of the most impor- 
tant results for the half-filling case within the present 
approach^, we extend the analysis of the mean-field solu- 
tions to the strongly correlated small-doping regime and 
discuss the evolution from weak to strong e-ph coupling, 
with particular attention to the way electronic properties 
are modified. The last Section is devoted to conclusions 
and acknowledgments. 



II. MEAN-FIELD ANALYSIS OF THE 
HUBBARD-HOLSTEIN MODEL 



The Hamiltonian of the Hubbard-Holstein model 
reads: 

h = -t ^2 c L c iv +uJ2 n ^ n a 

+ ujq a] dj + au) y)( a l + a ») n i W 

i i 

where t is a nearest-neighbor hopping term, c\ a (ci a ) is 
the creation (annihilation) operator for an electron with 
spin a at the site i while rii — rii a is the local den- 
sity operator. The operator a\ (a^) creates (destroys) a 
phonon at the site i with energy ujq and a parameterizes 
the coupling between electrons and local displacements. 
In what follows we will introduce the adiabaticity pa- 
rameter 7 = u>o/Cdt and the dimcnsionless parameter 
A = a 2 uJo/Cdt — a 2 j, which measures the strength of the 
e-ph interaction. Here Cd is a factor depending only on 
the dimensionality and the topology of the lattice defined 
by the expression of the non-interacting kinetic energy at 
half filling, \eq\ = Cat. 



As mentioned above, we treat the e-ph interaction by 
means of a variational Lang-Firsov transformation U = 
e s , with: 



S = a}^ [(rii) + fi (m - (rii))} (ai - a\). 



(2) 



fi are variational parameters which measure the cou- 
pling between phonon displacements and density fluc- 
tuations for any value of the adiabaticity parameter. 
We expect fi to be equal to one in the antiadiabatic 
regime, where the standard Lang-Firsov transformation 
applies, and equal to zero in the adiabatic regime where 
the Born-Oppenheimer approximation holds. Averag- 
ing H e ff = e~ s He s on the vacuum state of the trans- 
formed phonons we obtain a purely fermionic Hubbard- 
like model with phonon-dependent renormalizations of 
the kinetic and correlation terms and with a rescaling of 
the chemical potential, all of which are controlled by the 
variational parameters fi. 

The effective electronic model is then solved within the 
Kotliar-Ruckenstein slave-bosons technique^. For each 
site we introduce four bosons ei,di,pia-, representing the 
four possible states on site i (zero, two and one a electron 
state), and a fermionic operator Ci a which is connected to 
the original electron operator by the relation Cj CT = Zi a Cicr, 
where 
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The equivalence with the original Hilbert space is guar- 
anteed by the constraints 



4a c i* = pLp^ + 4 d i 

1 = ^2 PiaPio + d l d * ' 



(Vi,a) 
(Vi) 



which can be enforced introducing three Lagrange mul- 
tipliers a] 1 , A^/. The mean-field solution at a given 
value of density n in the paramagnetic homogeneous 
phase is obtained by taking the saddle-point value for 
the Bose fields ((e^) = eo, {di) — do, {pi a ) = Po) an d 
assuming translation invariance, so that fi = f and 
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Following closely Ref. |23|, we 



minimize the resulting variational energy with respect to 
the Lagrange multipliers and use the constraints to get 

E = -\e\qe- a2 f* +d 2 [U + 2a 2 io (f 2 -2/)] 

a 2 c (l - S) [1 - 5(1 - f ) 2 } (4) 

where e is the kinetic energy per site in the uncorrected 
case and n = 1 — 5; by introducing the standard notation 
x = eo + do 24 , on e can express the double occupancy 
as = (x 2 — <5) 2 /4x 2 and q = z 2 , i.e. the reduction 
of the kinetic energy due to the electronic correlation, 



as q 



{2x 2 



5 )/(! — 5 2 ) . In the absence of 
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e-ph coupling, when a = 0, Eq. J3J reduces to the well- 
known Gutzwiller energy for the pure Hubbard model 25 . 
On the other hand, taking the limit / = 1, one recovers 
the standard Lang-Firsov result, being the kinetic energy 
exponentially renormalized with the electron energy level 
and the Hubbard term shifted respectively by — a 2 ojo and 
— 2a 2 ojQ- 

In order to determine the mean-field solutions we min- 
imize (@J with respect to the remaining variational pa- 
rameters x 2 and /. One gets: 



e e 



2 fl 

-a f _ 



1-X 2 
1-5 2 

q\e\fe~ a2 f 2 =Wq(1-/) 



= [U + 2a 2 u; (f -2f)] 



(x 2 - 8f 
2x 2 



+ 8(1-8) 



- (5) 
• (6) 



For small doping 5 the kinetic energy can be expanded 
around the half- filling value \eo\ as M = ko|(l — a ^ 2 ) 
(with a depending on the specific shape of the uncor- 
related band), and the mean- field equation (JSJ can be 
conveniently rewritten as: 



(1 



-S 2 



1 
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(7) 



where u = [U/U c + A(/ 2 - 2/)/4] /e~ Q f , in which 
U c = 8|eq| is the Brinkmann-Rice critical value for the 
Mott transition in the absence of phonons. Eq. (7J coin- 
cides with the result for the Hubbard model once U/U c 
is replaced by u (cfr. Eq. (8) of Ref. l23h . This finding 
is easily interpreted in terms of renormalized interaction 
U e ff = U + 2(f 2 — 2f)a 2 u!o and renormalized kinetic 
energy e e ff = ee~° 2 f 2 , showing that the value of the 
parameter / determines to what extent the electronic 
properties are affected by phonons. Furthermore, being 
Eq. 10 a transcendental equation, it can be useful to 
separate the exponential from the algebraic dependence 
on /, obtaining: 



/ = 



2 1-aS 2 x 2 (2x 2 
1 + - 



2 r2 



7 1~8 2 x 4 + 8 2 (l-2x 2 ) 



(8) 



This form, as we shall see in the next Section, is more 
suitable for a graphical analysis. 

III. RESULTS 

We are now in the position to discuss the effect of e-ph 
interaction on the correlation-driven Mott transition at 
(5 = 0. The transition occurs for u = 1, which implies 
x = 0. In the small 7 adiabatic regime if we plug this 
condition into Eq. (JHJ), we find / = 7/4i£, while f = 1 
is recovered for large 7. Therefore, while in the antia- 
diabatic limit the e-ph coupling strongly reduces both 
the Coulomb repulsion and the kinetic term, in the adi- 
abatic regime, the small value of / makes the effect of 
e-ph smaller on both quantities. In this regime the expo- 
nential factor is closer to one, since the phonons are less 



effective in reducing the electron mobility, and, notice- 
ably, (U/e ) e ff = (tV £ o)(l — \l a2liJ o/U), in remarkable 
quantitative agreement with DMFT results^. For small 
7 the line marking the Mott transition in the A-U diagram 
is then given by A = (16/j)(U/U c — 1), and reproduces 
accurately previous diagrams obtained by DMFT 12 .. The 
increase of the critical U as a function of e-ph interaction 
results from the two renormalizations of U e ff and £ e //- 
Because both effects are essential to get the agreement 
with the DMFT result at 8 = 0, we expect them to be 
important also at finite 8. 

Let us now move away from half-filling, and consider 
the large U limit. Here we expect x 2 = \S\/( (being 
x 2 = at half- filling for u = U/U c > 1). Expanding Q 
to second order in 8 we obtain £ = y/1 — 1/u that again 
coincides with the expression for the Hubbard model 2 ^ 
with renormalized parameters depending on /. The value 
of / at small doping is then found by solving the self- 
consistent equation (derived from Eq. ©) 



/ = 



1 + B^e~~ s 

7 



(9) 



with B = 2u/(2u — 1). In the limit u ^> 1 (which is 
equivalent to u 3> 1 if A is not too large) the prefactor B 
goes to 1, thus simplifying the analysis of Eq. ©. 
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FIG. 1: Evolution with A of the numerical solution of the 
mean-field equations for u — 4 and various values of 7. 

It is useful to compare the behavior of Eq. © with the 
mean-field theory of the ferromagnetic transition. We 
will see that /, 7 and A play the roles of the magne- 
tization, temperature and magnetic field, respectively. 
The mean-field equation has the form / = h(f) and 
can be solved graphically just like the Curie- Weiss equa- 
tion. h(f) varies from to 1, so that / takes physically 
meaningful values. h(0) is different from zero for any 
finite value of 7 and h(l) < 1. Therefore h(f) inter- 
sects the straight line in at least one point. The function 
h(f) has an inflection point controlled by the value of 
the parameter 7. For small 7 the inflection point is at 
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small /, so that h(f) crosses / in three points, while for 
larger 7 only one intersection survives as the inflection 
point moves toward large /. 7 clearly plays a role similar 
to the temperature, with a critical value 7 C which sepa- 
rates the two regimes. The critical point (/ C ,7 C , A c ) may 
be evaluated analytically by imposing f c = h(f c ;X c ,Jc), 
1 = jf h Uc,K,lc), and = J»r/i(/ c ;A c ,7 c ), which 
in the limit B = 1 give the critical values f c = |, 
7c = 4e~l ~ 0.8925, and A c = ^ c ~ 3.0f23. When 
three solutions exist, we find two locally stable solutions 
/_ and /+ associated to weaker and stronger effective e- 
ph interaction. The two solutions correspond to negative 
and positive magnetization in the ferromagnetic termi- 
nology. Notice that the variational nature of our treat- 
ment of phonon degrees of freedom implies that only the 
lowest energy state is a valid result of our approach, even 
if the equations allow for more solutions. The parame- 
ter A (or more precisely the deviation from A c ) acts as 
a magnetic field in determining the existence of the two 
potential solution and which one is the ground state. For 
small (large) A only /-(/+) exists, and the energies of the 
two solutions cross in a first order line A c (7), ending in a 
critical point. To summarize, in the anti-adiabatic limit 
(large "temperature") regime, we have a single solution 
smoothly evolving with the "magnetic field" A. When 
A is large, lowering the "temperature" going from the 
anti-adiabatic to the adiabatic regime does not imply a 
dramatic change in the value of /. On the other hand 
for small A / changes from / w 1 to / w 7 moving from 
large to small 7. 

The above "Curie- Weiss" scenario is recovered also for 
finite yet large u = 4, by solving numerically the mean- 
field equations. In Fig^we plot / as a function of A for 
selected value of 7. As in the u ^> 1 limit, the evolution 
from weak to strong coupling is smooth for 7 > -f c ~ 0.9, 
while for 7 < j c the two solutions with small and large 
/ coexist. Again for large A / « 1 for any value of the 
adiabaticity parameter, while for small values of the e-ph 
coupling / tx 7 in the adiabatic regime and is / « 1 — 1/7 
when 7> 1. In Fig. |2Jwe report the first-order transition 
line between weak-coupling and strong-coupling solution, 
as well as the lines A c i and A C 2 marking the existence of 
the two solutions /+ and /_ . For large 7 the line of first- 
order transitions ends in a second-order end-point, where 
also A c i and A C 2 collapse. The main effect of a finite u 
compared to the infinite u limit is to move the critical 
point to slightly larger values of /, 7 and A. The strong- 
coupling solutions always show a full renormalization of 
U, while the kinetic energy is strongly suppressed moving 
from large to small values of 7. On the other hand, the 
weak-coupling solutions are characterized by a smooth 
renormalization of the correlation term that goes from its 
bare value in the adiabatic limit to the fully renormalized 
value U — 2a 2 uiQ as 7 — > 00; however the renormalized 
kinetic energy £ e ff, dominated by the exponential term 
exp(— Xf 2 /-f), is always rj 1. 

The phase diagram in Fig. [21 is not changed much as 
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FIG. 2: Phase diagram in the 7 — A space at u = 4. A c i and 
A C 2 mark the boundaries of the region in which respectively 
the /+ and /_ solutions disappear. A c (7) is the transition line 
where the two solutions have the same energy. 



function of doping. This result may be understood in the 
following way. The phase transition line is obtained by 
comparing the values of the energies for the solution for 
/_ and f + . As it can be inferred from Eq. (@J, both the 
kinetic energy and the polaron energy differ by a quantity 
of order 5 in the two solutions, so that the critical values 
of A depend on doping through subleading corrections in 
5 which are quite small at u = 4. The physics underlying 
the phase diagram of Fig. |2Jis that by increasing A a small 
number S of localized holes strongly coupled to phonons 
have a lower energy (~ — 5A|£o|, i.e., the energy of 5 
polaronic holes) than a bad metal that has a reduced 
kinetic energy, again of order S (e e // ~ — S\eq\) due to 
the large effective mass (to* ~ 1/6). 

We now consider the fate of the two solutions /_ and 
/+ as we approach the Mott transition by letting both 
(5^0 and U — > U c . In the adiabatic regime the mean 
field equations can be solved analytically and the pola- 
ronic solution only establishes at A po / = A(U/U e )(l — 
i/l — U c /U). At U ~ U c this gives A po ; = 4, which 
is reduced to 2 in the large U limit. We can translate 
these values of the dimensionless coupling in terms of 
the dimensional phonon-mediated attraction at the zero- 
frequency, Vph — 2A|eo|. We obtain = 8|eo| = U c 

close to the Mott transition and V^ 1 — 4|e | = U c /2 
in the large-[/ limit. We therefore find that, for U 
close to U c , the weak-coupling solution /_ is always 
the stable one even for quite large e-ph coupling up to 
Vph — U, thus showing that the Mott transition is robust 
with respect to e-ph interaction and polaron formation. 26 
Through the stabilization of the /_ solution the strong 
correlations close to the Mott transition also protect the 
quasiparticles from a strong mass renormalization due 
to phonon effects^. Considering the fluctuations around 
the mean-field solution one recovers the high-energy Hub- 
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bard bands^i. We expect that phonon effects wifl con- 
tribute significantly to themi4. 

As it is well known, the mean-field treatment may be 
incorrect in reproducing the order of a phase transition, 
and in this sense the polaronic transition discussed above 
may be turned to a crossover when fluctuations are in- 
cluded, as it happens in DMFT. On the other hand, the 
comparison with the DMFT results shows that the mean- 
field description captures the basic qualitative features of 
how the polaron formation is affected by the presence of 
the e-e interaction. 

We gain further information on the way the system 
evolves when the Mott insulator is approached by con- 
sidering the compressibility k = dn/d/i, or equivalently 
the Landau parameter Fq , related to k by the relation 
k = 27V*/ [1 + Fq], where N* is the quasiparticle density 
of states per spin at the Fermi level. To leading order in 
5 we obtain 



No/26 



4M 



2u- 1 



16|e |a 2 /(W) 



2u-l 
(10) 

where No is the free DOS (per spin) at the Fermi level. 
The second term in Eq. (|10|l represents the phononic 
contribution to the Landau parameter dressed by the e-e 
interaction. The first term describes the electronic con- 
tribution to Fq dressed by the e-ph correlation; once 
more we find the same expression of the Hubbard model 
with phonon-renormalized parameters^ 3 -. The attractive 
phononic contribution to Fq increases the compressibil- 
ity, and could eventually lead to a phase separatio n 13 ! 28 . 
This term is indeed only effective for the weak coupling 
/_ solutions. However the compressibility of /_ so- 
lutions diverges only for values of the parameters A, 7 
for which the /+ solution is energetically favored. This 
finding is closely reminiscent of the DMFT analysis of 
Ref. Il3l where the insulating behavior is favored by 
the e-ph interaction, but the divergence of the com- 
pressibility in the metal is prevented by the stabiliza- 
tion of the insulator. To verify the accuracy of our ap- 
proach, we compare our compressibility (|10|l to the pre- 
dictions of lar ge- N expansion^ and DMFT 13 . To com- 
pare with Ref. Ha . we take the large-U and small-7 limit; 
then H1(J|) acquires a simple random-phase approxima- 
tion form k(X) = 2AT*(0)/[1 + (F s ) e + (F s ) p/l ], where 
N*(0) = Nq/ 28 is the density of states per spin in the 
large-£7 limit in the absence of e-ph interaction and the 
phononic contribution to the symmetric Landau ampli- 
tude is simply given b y (F q ) p h — — 4A r *(0)a 2 w , in ex- 
cellent agreement with[lj| Relaxing the large-£7 limit we 
obtain 
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«W = 

k(0) 1 - 2a 2 uj Pn{0) 



(11) 



where k(0) is the compressibility in the absence of the 
e-ph interaction and j3 is a function of 7 that approaches 
unity from below in the adiabatic limit. Eq. is 
formally identical to Eq. (4) of Ref. IT3L where [3 is intro- 
duced as a fitting parameter in order to describe the de- 
viations from the random-phase-like expression for k(X) 
due to the finite (as opposed to infinite) value of U. 



IV. CONCLUSIONS 

In the present paper we have discussed an analytical 
approach to the subject of e-ph interaction in strongly 
correlated systems, based on slave bosons, treated at 
the mean-field level, and a variational Lang-Firsov trans- 
formation. The mean-field description shows that the 
Mott transition is robust with respect to e-ph interac- 
tions, yielding results in good quantitative agreement 
with previous DMFT studiesi*, and supplementing the 
DMFT scenario with some analytical results that make 
the physical picture more transparent. Furthermore we 
have shown that the nature of the evolution from weak 
to strong e-ph coupling is controlled by the adiabatic 
parameter 7. In the antiadiabatic regime, i.e., for 7 
larger than a critical value 7 C , the evolution is a smooth 
crossover, which is replaced, below j c , by a first-order 
transition where the "rotation" parameter / jumps from 
a small (/ ~ 7, almost no rotation) to a large value 
(/ ~ 1, fully developed polaronic regime). /, despite be- 
ing simply a variational parameter, measures the effec- 
tiveness of the e-ph coupling in affecting the electronic 
properties. 

In the adiabatic regime the phonon-mediated attrac- 
tion which screens the coulomb repulsion is found to 
be reduced proportionally to 7, in good agreement with 
DMFT analyses. In this regime the phonon-induced ef- 
fective mass enhancement is also reduced proportionally 
to the degree of adiabaticity. The first-order character of 
the transition from weak to strong coupling is most likely 
a spurious outcome of the mean-field theory. However the 
quantitative estimates of the boundaries between regions 
with different physical behavior are probably sound: at 
least they provide the correct order of magnitude in those 
cases where we have accurate estimates from alterna- 
tive approaches. To what extent the present approach 
can give reliable results far from the strongly correlated 
regime is an interesting open question but beyond the 
scope of this paper. 

We acknowledge discussions with S. Ciuchi and G. San- 
giovanni, and financial support from MIUR under Grant 
No.COFIN2003020239-006. 
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